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Abstract —Far-field Wireless Power Transfer (WPT) and Si¬ 
multaneous Wireless Information and Power Transfer (SWIPT) 
have attracted significant attention in the RF and communica¬ 
tion communities. Despite the rapid progress, the problem of 
waveform design to enhance the output DC power of wireless 
energy harvester has received limited attention so far. In this 
paper, we bridge communication and RF design and derive novel 
multisine waveforms for multi-antenna wireless power transfer. 
The waveforms are adaptive to the channel state information and 
result from a posynomial maximization problem that originates 
from the non-linearity of the energy harvester. They are shown 
through realistic simulations to provide significant gains (in terms 
of harvested DC power) over state-of-the-art waveforms under a 
fixed transmit power constraint. [] 

I. Introduction 

Wireless Power Transfer (WPT) via radio-frequency radia¬ 
tion is attracting more and more attention. RF radiation has 
indeed become a viable source for energy harvesting with clear 
applications in Wireless Sensor Networks (WSN) and Internet 
of Things (IoT) (3, ED. Interestingly, since WPT and wireless 
communication share the same RF medium to transfer power 
and information, the emerging field of Simultaneous Wire¬ 
less Information and Power Transfer (SWIPT) has recently 
attracted significant attention in academia 13, 0- 

The challenge with WPT is to find ways to increase the 
DC power level at the output of the energy harvester without 
increasing the transmit power. Most of the technical efforts 
have been devoted to the design of efficient energy harvesters 
(so-called rectenna) JT), El- Interestingly, the overall RF-to- 
DC conversion efficiency of the rectenna is not only a function 
of its design but also of its input waveform. However, the 
waveform design has received less attention 0-0. In 0, 
0, a multisine signal excitation is shown through analysis, 
simulations and measurements to enhance the DC power and 
RF-DC conversion efficiency over a single sinewave signal. In 
0, various input waveforms (OFDM, white noise, chaotic) 
are considered and experiments show that waveforms with 
high peak to average power ratio (PAPR) increase RF-to-DC 
conversion efficiency. Even though those papers provide some 
useful insights into the impact of waveform design onto WPT 
performance, there are many limitations in the WPT waveform 
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design literature: 1) there has not been any formal tool to 
optimize waveforms for WPT so far, 2) multipath fading (well 
known in wireless communications) has been ignored despite 
its tremendous impact on the received waveform at the input of 
the energy harvester, 3) the Channel State Information (CSI) 
is assumed unknown to the transmitter, 4) the transmitter is 
commonly equipped with a single antenna. 

In the rapidly expanding SWIPT community, the non¬ 
linearity of the energy harvester is not accurately modeled. It is 
indeed assumed that the harvested DC power is modeled as a 
conversion efficiency constant multiplied by the average power 
of the input signal to the energy harvester 0, 0- This was 
used so as to simplify the design of SWIPT but is unfortunately 
unrealistic 0-0. 

In this paper we tackle the important problem of waveform 
optimization for Multiple Input-Single Output (MISO) WPT. 
We consider WPT with a multisine waveform (due to its 
popularity in communication, e.g. OFDM) on each antenna 
and transmission over a multipath channel. We introduce a 
simple and tractable analytical model of the energy havester, 
accounting for its nonlinearity. Assuming perfect Channel 
State Information at the Transmitter (CSIT) can be attained 
(in a similar way as it is done in wireless communication 
systems), we formulate an optimization problem to adaptively 
change on each transmit antenna the multisine waveform as 
a function of the CSI so as to maximize the output DC 
current at the energy harvester. The globally optimal phases of 
multisine waveform weights are obtained in closed form while 
the locally optimal amplitudes are shown to result from a non- 
convex posynomial maximization problem subject to a power 
constraint, which can be formulated as a Reverse Geometric 
Programming and solved iteratively. The optimized waveforms 
adaptive to the CSI are shown to provide significant gains over 
state-of-the-art waveforms. 

Organization: Section HU introduces the WPT system model 
and section ED the analytical model of the energy harvester. 
Section IIVI tackles the waveform optimization and section [3 
evaluates the performance. Section [VI] concludes the work. 

Notations: Bold lower case and upper case letters stand for 
vectors and matrices respectively whereas a symbol not in bold 
font represents a scalar. Operator ||.||p refers to the Frobenius 
norm a matrix. £ {.} refers to the expectation operator. 


II. WPT System Model 

Consider the transmitter made of M antennas and N 
sinewaves whose multisine transmit signal at time t on transmit 
antenna m = 1 ,,M is given by 

JV-l 

x m(t) = ^ \ ^n,m COs)w n t + fin,m) (1) 

n= 0 

where s n/m and cf> n ,m refer to the amplitude and phase of the 
n th sinewave at frequency w n on transmit antenna m, respec¬ 
tively. We assume for simplicity that the frequencies are evenly 
spaced, i.e. w n = wo + nA w with A w the frequency spacing. 
The magnitudes and phases of the sinewaves can be collected 
into matrices S and $ such that the (n, m) entries of S and $ 
write as s n ^ m and respectively. The transmitter is subject 
to the power constraint J2m=i^{ \ x m \ 2 } = \ ||S||^ < P. 

The multi-antenna transmitted sinewaves propagate through 
a multipath channel, characterized by L paths whose delay, 
amplitude, phase and direction of departure (chosen with 
respect to the array axis) are respectively denoted as 77 , < 37 , ti 
and 0i, l = 1,..., L. We assume transmit antennas are closely 
located so that 77 , 07 and tl are the same for all transmit 
antennas (assumption of a narrowband balanced array) ED- 
The signal transmitted by antenna m and received at the single¬ 
antenna receiver after multipath propagation can be written as 

y (m \t ) (2) 

N—l / L—l \ 

= y ' Sn,m I y ' 07 COS {lV n [t Tl) ~\~ tl -p fi n ,m T" J 

n= 0 \ i—0 J 

where Arefers to the phase shift between the rn th 
transmit antenna and the first one. For simplicity, we as¬ 
sume that A n \ i = 0. For a Uniform Linear Array (ULA), 
A n ,m,l = 27r(m — lj^-cos ( 61 ) where d is the inter-element 
spacing, X n the wavelength of the n th sinewave. 

The quantity between the brackets in (0 can simply be 
rewritten as 

L —1 

COs{w n (t - Tl) + tl + (f> 

n,m H” 

1=0 

= COS {w n t ~\~ (3) 

where the amplitude A nyrn and the phase 'ijjn.m are such that 

A pj'4 , n,m _ A _ p j(pn.,mU (A\ 

- rL n,m c ' ° ll n,m V v 

with h n , m = = Ef^o 1 aiei(- w " T ‘ +A ^‘+M the 

frequency response of the channel of antenna m at w n . 

The total received signal comprises the sum over all transmit 
antennas, namely 

M M AT-1 

yit) = ^2 2/ (m) ( i ) =^2^2 Sn ’ mAn ’ m C0S ( W nt + V’n.m) 

m= 1 m= 1 n= 0 

N—l 

= ^2 Xn COS (W n t + S n ) (5) 

71—0 

where — \ e A pj'4 ), n,m 



Fig. 1. Antenna equivalent circuit. 
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Fig. 2. A single diode rectifier. 


III. Analytical Model of the Energy Harvester 
and DC Component 

In this section we derive a simple and tractable model of the 
energy harvester circuit and express the output DC current as 
a function of the waveform parameters. The model relies on 
several assumptions that are made to make the model tractable 
and therefore be able to optimize the waveforms. Performance 
evaluations will be conducted in Section using a more 
accurate circuit simulator. 

A. Antenna Equivalent Circuit 

Assume an energy harvester whose input impedance Ri n 
is connected to a receiving antenna as in Fig. Q] The signal 
y(t) impinging on the antenna has an average power P av = 
£{y(t) A ). Following J9j, the antenna is assumed lossless and 
modeled as an equivalent voltage source v s (t) in series with 
an impedance R an t = 50(1, as illustrated in Fig. |T] In the 
perfect matching case ( Ri n = Rant), the received power 
P av is completely transferred to the energy harvester’s input 
impedance such that P av = £{vi n (t) 2 }/Ri n where Vi n {t) 
is the input voltage to the energy harvester. Under perfect 
matching, Vi n (t) is half of v s (t) and both can be related to 
the received signal y(t) as 

tts (f) — 2 y(t) \j Rant, Vin {t) — y)t ) \J Rant • ( 6 ) 

B. Rectifier and Diode Non-Linearity 

A rectifier is always made of a non-linear device (e.g. diode) 
followed by a low pass filter (LPF) with load 13, 0, j6j. A 
simplified rectifier circuit is illustrated in Fig. [3 We assume 
that its input impedance has been perfectly matched to the 
antenna impedance. 

The current id(t) flowing through an ideal diode (neglecting 
its series resistance) relates to the voltage drop across the diode 

v d (0 

Vd(t) = v in {t) - v out (t) as i d (t) = i s (e^F - l) where i s is 
























the reverse bias saturation current, vt is the thermal voltage, 
n is the ideality factor (assumed equal to 1 for simplicity). 
In order to express the non-linearity of the diode, we take a 
Taylor expansion of the exponential function around a fixed 
operating voltage drop Vd = a such that the diode current can 
be equivalently written as 

OO OO 

id{t) = ^2 ki{v d {t)-a) 1 = ^2 ki(vin(t)-v out (t)-a) 1 , ( 7 ) 

2=0 2=0 

where k 0 = i s (e ^ - l) and k t = is ^Zl) 1 ’ i = 

As such, it is not easy to infer from i[7Ji the exact dependencies 
of the diode current on the waveform parameters since both 
Vin(t) and v ou t(t) will depend and fluctuate over time as a 
function of the waveform. Nevertheless, assuming a steady- 
state response, an ideal rectifier would deliver a constant 
output voltage v ou t that would track the largest peaks of the 
input voltage Vi n {t) j9]. As a consequence, the output current 
delivered to the load i out would also be constant. Denoting 
the magnitude of the peaks of Vi n (t ) as {£„, v ou t = v ln ■ In 
this ideal rectifier, since v ou t is a constant (we drop the time 
dependency), a suitable choice of the operating voltage drop 
a would be a = £[v in (t) - v out } = —v out = —v in since 
£ {Vi n (t)} = \JR a nt£ {j/(f)} = 0. Under such assumptions, 
<0 can simply be written as id(t) = ^i^antV^Y > 

which now makes the dependency between the diode current 
id(t), the received waveform y(t) and therefore the transmitted 
waveforms (x m (t)} much more explicit. Moreover, with such 
a choice of a = —v ou t, the Taylor expansion of the exponential 
function could be truncated without compromising too much 
on the accuracy. 

C. Output DC Current 

The problem at hand will be the design of {x m (f)} such that 
the output DC current is maximized. Under the ideal rectifier 
assumption, the current delivered to the load in a steady-state 
response is constant and given by i out = £ {id(t)}. 

In order to make the optimization tractable but keep the 
fundamental non-linear behaviour of the diode, we truncate 
the Taylor expansion to the fourth order, as argued in © to 
be sufficient to demonstrate the rectification operation. The 
output current therefore approximates as 

4 

iout = £ {id(t)} « ^2 kiR % ilt£ {y(t) 1 } . ( 8 ) 

i=0 

Applying <(5j to ([Sj and taking the expectation over time, we 
get an approximation of the DC component of the current 
at the output of the rectifier (and the low-pass filter) with a 
multisine excitation over a multipath channel as i out = fco + 
zdc where 

ZDC = k 2 Rant£ {y(t) 2 } + k 4 R 2 ant £ {y{t) A } (9) 

is detailed in (ITOt (at the top of next page). There is no third 
order term since £ {j/(t) 3 } = 0. Expression (j9j was used 


in |6) for single-antenna WPT to show the suitability of in- 
phase multisine excitation with uniform power allocation for 
scenarios without multipath and in the absence of CSIT. 

It is worth contrasting © with the model commonly used in 
the SWIPT literature. In 0, 0, the total harvested power is 
defined as E = (£ {y(f) 2 } where £ is a constant that refers to 
the RF-to-DC conversion efficiency and is commonly chosen 
as 1 for simplicity. Comparing with ©. we note that E only 
accounts for the second order term in the Taylor expansion and 
ignores the non-linearity of the rectifier (modeled by the fourth 
order term). The maximization of E subject to the transmit 
power constraint would lead to a single-sinewave transmission 
strategy where the power is allocated to the strongest sinewave, 
i.e. the one corresponding to the strongest channel lUOl . This 
contrasts with RF experiments that highlight the usefulness of 
allocating power to multiple sinewaves 0-0- If non-linearity 
is accounted for, the presence of the fourth order term in <[9]» 
suggests that transmitting over a single sinewave is in general 
suboptimal. 

IV. Waveform Optimization 

Assuming the rectifier characteristics ko,k 2 ,k 4 and the 
channel state information (in the form of frequency response 
h n ,m) are known to the transmitter, we now aim at finding 
the optimal set of amplitudes and phases (across antennas and 
frequencies) S, $ that maximize i out , i.e. 

max (11) 

S,*i> 

subject to — ||S||^ < P. (12) 

From the previous section, we however note that fc 0 , k 2 ■ /m 
are functions of a = —v ou t which is affected by the choice of 
S, 3>. Similarly, S, <f> are functions of ko,k 2 , k/\. This suggests 
that ko : k 2 , &4 and S, $ should be iteratively computed by 
fixing k 0 ,k 2 ,k 4 when S, $ are optimized and inversely. We 
here assume that ko,k 2 ,k 4 have been computed for a given 
a and we aim at finding the optimal S, <I>. For a fixed a, kg 
is fixed and is not affected by the choice of S, $. Hence the 
problem (ITU can equivalently be re-written as 

max z D c{ S,$) (13) 

S.<I> 

subject to i ||S||^ < P. (14) 

Interestingly, the optimal phases can be obtained first in closed 
form. Given the optimal phases, the optimal amplitudes can 
be computed numerically. 

A. Phase Optimization 

To maximize z/jc(S. d>), we should aim at guaranteeing 
that all cos(.) are equal to 1 in ( I Kil l. This can be satisfied by 
choosing tl> nym = 0 Vn, m (and therefore S n = 0 Vn), which 
implies from © to choose the optimal sinewave phases as 
= — i>n,m ■ is obtained by collecting m Vn, m 

into matrix. 




( 10 ) 


2 DC (S,^) = y-Rant 


E E Sn,mo Sn,mi -^n,mo -An,mi COS (lpn,mQ ' 1 pn,mi ) 

n =0 mo,mi 


3^4 2 

0 ^ant 


E 

E 

' 3 

IE! Sn i’ m i A nj,mj 

COS (071q ,777-0 H - 077. 1 ,777, 1 '077-2 > m 2 077-3,777.3 ) 

n 0 > n l > n 2 > n 3 

_ 710+72^=712 + ^3 
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Li=o 
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B. Amplitude Optimization 


{ 7 /c}, problem (1 1 7li-( l 1 9b is now replaced by the standard GP 


With the optimal phases, i l> n ,m — 0 and X n = 

such that Adc(S,$*) is simply obtained 
from m by replacing all cos(.) functions by 1. 

Recall from ED that a monomial is defined as the function 
g : R^ + — > R : g(x) = cx^x^ 2 ■ ■ -x 0 ^ where c > 0 and 
£ 1. A sum of K monomials is called a posynomial 
and can be written as /(x) = XE =1 <?fc(x) with c/fc(x) = 
Ckx\ lk X 2 2k ■ ■ ■ x°^ k where Ck > 0. As we can see from (ITOt . 
zdc(S, $*) is a posynomial. 

The optimization problem becomes 

max z DC ( S,$*) (15) 

subject to i ||S||^ < P. (16) 

It therefore consists in maximizing a posynomial subject to 
a power constraint (which itself is written as a posynomial). 
Unfortunately this problem is not a standard Geometric Pro¬ 
gram (GP) but it can be transformed to an equivalent problem 
by introducing an auxiliary variable to 

min 1/to (17) 

S,£o 

subject to i ||S||^ < P, (18) 

ZDc(S,^)/t 0 >l- (19) 

This is known as a Reverse Geometric Program due to 
the minimization of a posynomial subject to upper and 
lower bounds inequality constraints UTI . ill2l . Note that 
ZDc(S,&*)/to > 1 is equivalent to fo/z.Dc(S, <F*) < 1. 
However 1 /tdc(S, d?*) is not a posynomial, therefore pre¬ 
venting the use of standard GP tools. The idea is to lower 
bound z/x;/(S. <!>*) by a monomial z/x::(S), i.e. upper bound 
1 /^dc(S,^*) by the monomial 1 /tdc(S) (since the inverse 
of a monomial is still a monomial) G2, Q3).Let {g k ( s,$*)} 
be the monomial terms in the posynomial £dc(S, T**) = 
^2k=i •!**)■ The choice of the lower bound relies on the 

fact that an arithmetic mean is greater or equal to the geometric 
mean. Hence, 2 dc(S,$*) > Tlfeli (flfc(S, &*)/'Yk) lk = 
zdc( S), where 7 fc > 0 and Yjk= l = 1- Since 

1/z dc (S,**)<1/zdc(S), (20) 

we can replace (in a conservative way) in¬ 
equality t 0 /z DC { S,$*) < 1 by t 0 /z DC {S) = 
f 0 nf=i (5fe(S,$*)/7fc) _7it < 1- For a given choice of 


min 

s,t 0 


subject to 


1/io 

\ II s IIf < 


n (M 

*:=i v ’ 


~ik 

< i, 


( 21 ) 

( 22 ) 

(23) 


that can be solved using existing software, e.g. CVX Ga¬ 
it is important to note that the tightness of the upper bound 
(l20l > heavily depends on the choice of { 7 /.}. Following lITTIl . 
m, an iterative procedure can be used where at each iteration 
the standard GP (l2TI) - (l23l) is solved for an updated set of 
{ 7 k}- Assuming a feasible set of magnitude S'7 at iteration 
i, compute Tfe(sW) = ^*)/z DC (S {i \ $*) VAs. Then 

solve problem (l2Tl i-(l23l> to obtain S( I+1 E Repeat the iterations 
till convergence. 

As noted in mi, since the original problem is nonconvex, 
the final solution is not guaranteed to be the global optimum. 


V. Simulation Results 

We now evaluate the performance gain of the optimized 
adaptive waveform versus two baselines: a non-adaptive wave¬ 
form not relying on CSIT and an adaptive waveform relying on 
CSIT but not requiring the optimization of Section [IV] A suit¬ 
able choice of non-adaptive waveform for single antenna WPT 
was shown in 0 , 0 to exhibit high PAPR. It was suggested 
to choose an in-phase multisine excitation with uniform power 
allocation. We therefore choose the non-adaptive baseline 
waveform as (f> n ,m = 0 and s n)TTl = l/y/NM \/n,m. The 
adaptive baseline waveform is chosen as a matched filter (MF) 
allocating power to all sinewaves but proportionally to the 
channel strength, i.e. (+,, m = and s n>m = cA n ^ m with 

c a constant to guarantee the power normalization. Hence the 
difference between the optimized waveform and the one based 
on MF lies in a different choice of amplitudes, hi = 0.0034 
and &4 = 0.3829 have been computed for an operating point 
a = 0 and used as such to design the optimized waveform. 

We assume a WiFi-like environment at a center frequency 
of 5.18GHz with a 36dBm EIRP, 2dBi receive antenna gain 
and 50dB path loss in a large open space/office environment 
with a NLOS channel power delay profile obtained from 
model B m. Taps are modeled as i.i.d. circularly symmetric 
complex Gaussian random variables. This leads to an average 
received power of about -12dBm. The frequency gap is fixed 
as A w = 2nAf with Af = B/N with B = 20 MHz and the 
N sinewaves are centered around 5.18GHz. 













Fig. 3. Average zjjc as a function of (IV, M). 


Fig. 5. Average DC power delivered to the load as a function of (IV, M). 
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gains (in terms of harvested DC power) over state-of-the-art 
waveforms under a fixed transmit power constraint. The results 
are expected to trigger significant interests in the RF/WPT 
community as well as in the communication theory community 
involved with wireless EH communication and SWIPT. 
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Fig. 4. Rectenna with a single diode and a L-matching network. 
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